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THE MEASUREMENT OF SPEED IN MOTOR PERFORMANCE* 


C. H. McCLoy 


STATE UNIVERSITY OF IOWA 


When the centroid method of factor analysis was applied to two 
sets of data on athletic performances, three significant factors 
emerged: strength, velocity, and dead weight. Scores on this speed 
factor were predicted by the multiple regression technique, the factor 
loadings on the speed factor being used as the criterion correlations, 
and these predicted scores were correlated with each of the other 
variables. When the original tables, augmented by the new speed 
variable, were refactored, the computed speed factor fell on the speed 
axis as a primary trait. It is thus shown that it is possible to isolate 
and measure a factor which appears in variables under consideration 
only as a compound. 


In the analysis of many of the simpler motor skills used in physi- 
cal education and athletics, such as in the runs, the jumps, and the 
throws, three orthogonal factors usually emerge, which are readily 
identifiable as (1) strength, or the force of contraction of muscular 
tissue, (2) speed of movement, and (3) dead weight. For the pur- 
poses of preliminary orientation, it may be well to discuss these fac- 
tors briefly before considering the results of the study presented in 
this paper. 

In research in this type of motor activities, it is quite practicable 
to insert specific variables of relatively known composition. For ex- 
ample, there are available numerous tests of strength which involve 
no significant speed component but which involves simply slow, force- 
ful contractions. They are so simple that a few rehearsals suffice to 
enable the subject to master the skill; hence, educability in this spe- 
cific performance is not an appreciable problem. If a number of such 
tests are used in a factor study, they invariably cluster together in 
factorial space. (For example, see Table II, variables 6 to 9.) 

Speed, on the other hand, is not so readily definable. It is speed 
of movement and is orthogonal to strength. It must not be confused 
with power. In such events as a high jump, the velocity of projection 
upward from the ground is a resultant of .1e power expended. This 
power (P) is the product of force (F’) and velocity (V), (P = FV). 
We are, therefore, used to associating our concept of speed with this 

* Much of the expense of the statistical computations connected with this 


study was financed by a grant from the Carnegie Foundation for the Advance- 
ment of Teaching. 
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resultant of the expenditure of power. This may be made clearer with 
an illustration. Suppose we were to choose 1000 men of the same age, 
height, weight, and strength. They would not all respond, in motor 
athletic performance, with the same speed. Residuals from predicted 
velocity in actual samples (predicted in adults from height, weight, 
and strength) distribute themselves approximately normally. It is 
these differences in speed that are measured by this speed factor; and 
they may be due to differences of viscosity in the sarcoplasm of the 
muscular tissue, or they may be due to differences in leverage in the 
attachments of the muscles. 

For example, in one individual, the distance from the olecranon 
process of the elbow to the center of the joint with the humerus is 1.7 
cms. The distance from this point to the distal end of the ulna is 26.3 
cms. This relationship of the power arm to the work arm differs sig- 
nificantly in different men, resulting in greater speed at the distal 
end of the forearm with a given speed of movement at the end of the 
power arm. Whatever the cause, the result is measurable through a 
factor analysis. 

The third factor is dead weight. From an athletic point of view, 
weight is divisible into two parts. “Live weight” is that weight which 
is muscular and which pays its way. It is indirectly measurable 
through the strength factor. This, in the adult male, makes up about 
43 per cent of the total. “Dead weight” is the remainder—the freight 
that is carried by the muscles. This dead weight is positively corre- 
lated with weight-throwing events such as the shot put, for it aids 
through its inertia in resisting the reaction of the throw. It is fre- 
quently negatively correlated with the runs and jumps, for here it is 
just extra baggage. 

Since these three factors are so universally found in motor events, 
variables calling for them can easily be included in any appropriate 
battery of motor tests in order to isolate and identify these factors 
from among all of the factors that may emerge in the more compli- 
cated motor events. 

Since the track and field or power type of events falls on a plane 
that is oblique to that of the strength factor, it is plain that these 
events are compounds of strength and speed. Hence, we wished to 
isolate this speed factor for experimental purposes. 

In the present study, two sets of data were used. The first con- 
sisted of athletic performances of 100 college men for the following 
events: 


4-pound shot 
8-pound shot 
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12-pound shot 

16-pound shot 

24-pound shot 

Putting a baseball (in the same manner as for the shot) 

The unweighted strength index (total strength) 

The “athletic” or weighted strength index (has some 
dead weight) 

The Sargent Jump (springing straight upward) 

Weight 


These data were collected by Mr. James W. Coleman* of the Uni- 
versity of Nevada. The ages of the group ranged from 18 to 24. 

The second set of data consisted of athletic performances of 163 
junior-high-school girls in the following events: 


A short dash 

The standing broad jump 

Basketball throw for distance 
3-pound shot put 

12-pound shot put 

Back strength 

Leg strength 

Shoulder strength 

Grip strength 

Total strength 

The Sargent Jump 

Weight 


These data were collected by Miss Jane E. Harris} of Sioux City, 
Iowa. In both studies, the shot-putting events were done from a 
stand, without the usual preliminary hop. 

These data were factored and rotated by the “two-at-a-time” 
method, and three significant factors emerged, which were readily 
identifiable as the three mentioned above—strength, velocity, and dead 
weight. The rotated factor loadings are given in Tables I and II in 
roman type at the upper left. 

We then used multiple regression equations to predict the speed 
factor, using the factor loadings on the speed factor for the correla- 
tions with the criterion, and the intercorrelations between the vari- 
ables for the other correlations (see Tables III and IV). 

The resultant multiple regression equations for the speed factors 
were as follows (the results are in terms of 7-scores) : 


* Coleman, James W. The differential measurement of the speed factor in 
large muscle activities. Research Quarterly, 1937, 3, 123-130. 

+ Harris, Jane E. The differential measurement of force and velocity for 
junior high school girls. Research Quarterly, 1987, 4, 114-121. 
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For the men: 


Speed factor = 1.573 (4-pound shot in ft.) + .308 (Sar- 
gent Jump in cms.) — .0185 (Strength Index in 
Ibs.) — 2.512; 
R= .887. 


For the girls: 


Speed factor = .824 (Standing broad jump in in.) + 
.893 (3-pound shot in ft.) — .0505 (Strength In- 
dex in lbs.) — 7.00; 

R = .909. 


It will be noted that in the formulas above, the Strength Index,* which 
has approximately a zero correlation with the speed factor and with 
the dead weight factor, was used to partial out the strength in the 
power events used. 

Since a factor has no absolute mean or standard deviation,} we 
used 50 for the mean and 10 for the standard deviation, and pre- 
dicted the factor in the form of a T-score.t The predicted speed score 
was found for each subject and correlated with each of the other vari- 
ables, and the material was refactored, yielding the centroid factor 
matrices in Tables V and VI. 

The resulting rotated factor loadings are found in Tables I and 
II in italics at the lower right. In determining the number of factors 
to which to carry the results, two criteria were used. The fourth fac- 
tor was computed and plotted against each of the others. The clusters 
seemed very slender along the axis and were entirely meaningless to 
the experimenter. Second, the third factor residuals were surveyed, 
and in view of the number of cases, seemed not to justify further com- 
putations. The distributions of the absolute values of the residuals 
were as follows: 


* The “Strength Index” is a total strength score made up of the sum of a 
number of strength test variables. See McCloy, C. H. Tests and Measurements in 
Health and Physical Education. New York: F. S. Crofts, 1939, Chapter 4. 

+ All of the other variables used have absolute values, regardless of the 
group studied. Thus a 10-year old boy and a 20-year old man, each of whom runs 
100 yards in 12 seconds, achieve the same record. The scores for a speed factor, 
however, computed for each of the two age groups, would probably be different 
since they would be in terms of standard deviations from means of different ex- 
perimental groups. 

t We prefer T-scores to z-scores (or standard scores) because we find them 
to be much more meaningful to students, who frequently seem to have a high 
natural resistance against minus signs. 
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Men Girls 

16 2 

14 1 

12 

10 1 
.08 1 4 
.06 6 6 
04 7 17 
.02 10 27 
.00 28 23 


These are the third factor residuals of the second process of factor- 
ing. The results are interesting. It will be noted that (Tables I and 
II): 


TABLE I 


Rotated Factorial Matrix 
College Men 


I II Ill h 

Strength Speed Dead weight 
factor factor factor 
.480 -770 077 910 

1. 4-lb. shot put 395 850 £27. 954 
489 641 476 936 

2. 8-lb. shot put 32 589 A455 915 
500 611 458 913 

8. 12-lb. shot put 477 501 611 923 
461 703 468 .962 

4, 16-lb. shot put ALS 583 633 969 
517 521 491 883 

5. 241b. shot put 494 420 651 919 
823 556 186 669 

6. Baseball put .262 551 -250 660 
.936 .000 .000 .936 

7. Strength index 949 030 004 950 
.890 —.021 241 .922 

8. Athletic strength 

index 978 228 —.136 1.013 

.288 577 —.326 722 

9. Sargent Jump 064 684 —.063 -690 
653 .010 431 -782 

10. Weight 703 —.022 317 771 

11. Computed Speed 

factor 082 930 —.003 934 


Figures in roman type, upper left of each space, are for the first factorial 
matrix. Figures in italics, lower right of each space, are for the second factorial 
matrix. The computed speed factor was computed from the first factorial matrix 
and hence was not included in it. 
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(1) The computed speed factor falls on the speed axis and has 
a correlation of .0 with the strength and the dead weight factors. 

(2) Its vector is approximately the same as the multiple corre- 
lation with the factor, as would be expected theoretically. The minor 
departures from the theoretical values are readily accounted for by 
the differences in the assumed values in the diagonals in correlation 
matrices of this size. The Strength Index measures the strength fac- 


tor in almost pure form. 


TABLE II 


Rotated Factorial Matrix 
Junior High School Girls 





I II Ill h 
Strength Speed Dead weight 
factor factor factor 
279 .642 —.150 -716 
1. Dash 347 578 —.101 677 
O27 775 —.265 .882 
2. Standing broad jump 821 854 —.256 947 
.557 613 .190 .850 
3. Basketball throw 
for distance .610 543 250 858 
555 .602 139 .830 
4. 3-lb. shot put 578 607 171 855 
.671 .421 352 .867 
5. 12-lb. shot put 700 400 347 877 
823 .051 .188 .846 
6. Back strength 797 031 119 806 
.862 —.056 —.145 .876 
7. Leg strength 873 —.208 —.033 896 
.620 .808 257 738 
8. Shoulder strength 627 SOL 187 722 
-760 .167 142 -791 
9. Grip strength 767 198 —.147 806 
1.019 .000 .000 1.019 
10. Total strength G95 —.020 —.051 996 
.330 -729 —.219 .829 
11. Sargent Jump B45 739 —.170 833 
519 .000 .524 -738 
12. Weight 629 —.152 823 728 
13. Computed Speed 
factor —.004 877 —.005 877 


Figures in roman type, upper left of each space, are for the first factorial 
matrix. Figures in italics, lower right of each space, are for the second factorial 
matrix. The computed speed factor was computed from the first factorial matrix 


and hence was not included in it. 
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The major contributions of the study seem to be the follow- 
ing: 

(1) Itis found possible to isolate and measure in relatively 
pure form, and with a high degree of accuracy, a factor which 
appears in the variables studied only as a compound (in oblique 
relationship). This is done by computing a multiple correlation, 
with 50 assumed as the mean and 10 as the standard deviation 
for the factor. 

(2) When replaced in the correlation matrix, this new var- 
iable (the speed factor) falls on the axis as a primary trait. 

(3) The variable may be used, practically, to compute, in- 
directly, the potential movement speeds of the body parts for this 
type of motor activity.* 

(4) This technique may be used to isolate other factors in 
this or in other factorial matrices. The factor loadings are used 
as the correlations of the independent variables with the factor, 
and the variables selected should consist of one or more that are 
heavily weighted with the factor to be measured, and one or more 
with zero loadings with that factor, but with relatively high cor- 
relations with the other variables used and with the other factors, 
in order that the influence of these other factors may be par- 
tialed out. 


* It has been so used by one coach of track and field athletics, who found that 
all boys showing a large amount of this factor (T-scores 65 and up) made excel- 
lent track athletes when trained, regardless of their previous experience. 
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TABLE V 


Centroid Factor Matrix 
(Second factoring) 


College Men 


I II 
4-pound shot .882 .276 
8-pound shot .903 .080 
12-pound shot .865 .109 
16-pound shot -902 .200 
24-pound shot .840 .048 
Baseball put .623 .210 
Strength index .636 —.653 
Athletic strength index -730 —.569 
Sargent Jump .488 .295 
Weight 591 —.467 
Speed factor .665 -508 

TABLE VI 


Centroid Factorial Matrix 
(Second factoring) 


Junior High School Girls 


Dash 584 252 
Standing broad jump .695 -450 
Basketball throw, distance .824 .183 
8-pound shot 823 .231 
12-pound shot .834 .046 
Back strength -703 —.363 
Leg strength .636 —.615 
Shoulder strength -708 —.035 
Grip strength -741 —.284 
Total strength .828 —.554 
Sargent Jump .663 369 
Weight .483 —.368 


Speed factor 454 -708 


III 
—.252 
121 
316 
302 
391 
—.034 
—.103 
—.294 
—.018 
306 
—.433 


—.231 
—.459 
131 
036 
-269 
158 
147 
138 
—.144 
.024 
—.343 
—.394 
—.247 
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WEBER LAWS, THE WEBER LAW, AND 
PSYCHOPHYSICAL ANALYSIS 


ALSTON S. HOUSEHOLDER AND GALE YOUNG 
THE UNIVERSITY OF CHICAGO 


In failing to define the units in which the stimulus is to be 
measured, the Weber law might seem to make no definite assertion, 
and indeed, it is shown that any single empirical function, supposed 
to relate a given stimulus intensity with that intensity which is just 
noticeably greater, can be put into the Weber form by a suitable 
change of scale in which the stimulus intensity is to be measured. 
Nevertheless, it turns out that if different individuals have different 
Weber functions, when the intensities are measured on a given scale, 
then it is by no means always possible to transform the scale so that 
all of the functions can take on the Weber form. Some necessary con- 
ditions are given for the possibility of such a transformation when 
there is at hand a finite number of functions, and when the func- 
tions depend upon a single parameter the necessary and sufficient 
condition is easily derived. The same discussion leads to a general- 
ization of Thurstone’s psychophysical scale and shows that such a 
scale is always possible. 


In discussing the discrimination of intensities there are two em- 
pirical questions that are not always clearly distinguished in the lit- 
erature. On the one hand, to a graded stimulus we may have a graded 
response, in the sense that if a stimulus of intensity S is applied a re- 
action of intensity R results, and the result of increasing the inten- 
sity of the stimulus is to increase the intensity of the response with- 
out otherwise altering it. Thus experimentally R is obtainable as a 
function of S which is presumably monotonically increasing and con- 
tinuous. Then one could choose some fixed quantity A R and ask for 
that value of A S associated with each S such that the responses 
elicited by S and by S+ 4S differ in the amount 4 R. Then by keep- 
ing A R fixed but varying S one obtains A S as some function of S. A 
function of this type could be taken as giving the discriminating pow- 
er of the organism, and frequently is so taken when studies are made 
of tropisms in lower organisms (3). 

But in the study of the higher vertebrates, at least, a somewhat 
different problem is usually set. Here stimuli of different intensities 
may lead to responses that differ qualitatively rather than quantita- 
tively—a rat can be trained to jump to the right on seeing a light of 
given brightness, and to jump to the left when the light is of a differ- 
ent brightness. Then the intensities are said to be “discriminated”, 
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or the difference to be “noticeable” by the organism, provided the 
responses called forth are qualitatively different. Then one can ask 
this question: If a stimulus of given intensity S calls forth a given 
response, by what amount A S must the intensity be increased in or- 
der to elicit a qualitatively different response? (It is understood that 
that A S is to be taken for which the response-difference occurs in a 
suitable, preassigned percentage of trials.) Again A S is empirically 
obtainable as a function of S , but a priori no relation can be predicted 
between this function and the one defined above. In the literature the 
two functions are not always clearly distinguished, although the very 
possibility of obtaining qualitatively different responses to quantita- 
tively graded stimuli implies the presence in the nervous system of 
some kind of a “straining” mechanism, not required for merely quan- 
titatively graded responses, which has the property of diverting exci- 
tations of different intensities along different nervous pathways (see, 
e.g. 5, 6, 8, and 9). Perhaps the reason for the confusion can be traced 
back to Fechner, who supposed that there was some intermediate 
psychic process of measure RF operating in discrimination of the sec- 
ond type, which is such that 4 R must exceed a certain fixed value 
before such discrimination can be made, i.e. before the overt response 
can undergo a qualitative change. 

Whatever the reason may be, the two functions are distinct, and 
we are here interested in the function of the second type. Thus, for 
any sense modality, and for any given choice of physical units in which 
to measure the intensity S of the stimulus, there is an empirically de- 
terminable function A S of S such that a change in S in the amount of 
A S brings forth a qualitative change in the reactions on the part of 
the organism being investigated. Otherwise stated, for any intensity 
S, there is a minimal stronger stimulus S. whose value is a function 
of S, 

S, = f (S,) ’ (1) 


for which the response changes qualitatively in a given percentage of 
the trials. Any such function f(S) so determined we shall call a Web- 
er function, and if, in particular, f(S) is linear and homogeneous: 


f(S) =as, (2) 


where a is a constant, then the function will be said to have the Weber 
form, or to be, simply, linear. Note that this constant a exceeds by 
unity the ordinary “Weber constant” A S/S. 

Weber’s law asserts, as is well known, that the Weber function 
has the Weber form. But Weber’s law fails to define the units in 
which the stimulus is to be measured, whereas obviously, if the Weber 
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function has the Weber form with one choice of units, it will not nec- 
essarily have this form when the units are changed. Thus suppose 
that for some choice of units S the Weber law is verified empirically. 


Then if one should introduce new units S defined by 


S=logS, S=e 
we have ¥ 
S.=a+5S., a = loga, 


so that the associated Weber function has no longer the Weber form. 
Thus with any arbitrary choice of units for the measurement of the 
intensity of the stimulus of a given kind, there will be associated a 
particular Weber function whose form is empirically determinable. 
If the function does not have the Weber form, it may be that by a 
new choice of units the associated Weber function would turn out to 
have this form. Hence we may ask the question: Given any Weber 
function, supposed to have been determined empirically after fixing 
upon a definite scale of units, what change of scale can be employed 
in order that the Weber function associated with the new scale units 
shall have the Weber form? If it should turn out that for an arbitrary 
function f(.S), supposed to be a Weber function, it is always possible 
to introduce a scale transformation which will give the transformed 
Weber function the Weber form, then the Weber law is not a law at 
all, and is of no value unless it be taken as a definition of the scale 
units or a class of scale units, viz. that scale in which the associated 
Weber function has the Weber form. But if it should turn out that 
functions f(S) exist of the type that might a priori be Weber func- 
tions but such that they cannot be transformed to the Weber form by 
any choice of scale units, then the Weber law can be made definite if 
it be taken as affirming that in fact no empirical Weber function will 
ever be one of these functions f(S). This statement could be verified 
or refuted empirically, at least if we may be permitted to disregard 
errors of observation. 

This is the problem for consideration, and we now define it more 
precisely. Let x or y, instead of S, represent the measure of the 
stimulus intensity in any scale. Let 


x =f (x) (3) 


or 
y=g9(y) (4) 
denote the measure of that stimulus which is judged just greater than 


x or y in a suitable, preassigned percentage of the trials. Evidently 
always 
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f(x) >2, g(y) >y, (5) 


at least if this “discriminal” percentage is taken large enough, as we 
shall suppose for the present. We suppose, moreover, that these func- 
tions are properly monotonically increasing, continuous functions of 
their respective arguments. Finally it is assumed that the “intensity” 
is an essentially positive quantity, so that arguments and functions 
assume only positive values, and that zero intensity is associated with 
the absence of the stimulus. 

Evidently the Weber functions f(2) or g(y) can be defined over 
a finite range of values of x or y only: too heavy a weight cannot be 
lifted, too bright a light is injurious, etc. Hence on any scale there 
will be a maximal x or a maximal y for which f(x) or g(y) is defined. 
It may be that at this maximal value the corresponding inequality 
(5) must be replaced by an equality. 

There will also be a minimal x or y, different from zero for which 
f(x) or g(y) is defined, because of the presence of a finite threshold. 
One could define f(0) or g(0) as the threshold value, which would be 
equivalent to including the absolute threshold among the differential 
thresholds, but this we shall not do. Hence f(x), or g(y), is defined 
for only those values of x or y that are supraliminal. 

Now suppose a given scale is chosen, the x-scale, and that the 
associated Weber function f(x) is determined empirically. Let g(y) 
be a function of prescribed form, arbitrary except for the restrictions 
imposed above, and let us inquire whether there exists a scale trans- 
formation 


Y= $¢(2) (6) 


which is such that the function f(x) transforms into g(y) by this 
scale transformation. We suppose that ¢(x) is essentially positive, 
and is properly monotonically increasing and continuous in x, and 
that it is defined over the range for which f(x) is defined. Naturally 
its values, as x varies within this range, must all be included in the 
range of values of y over which g(y) is defined. 

If the function ¢ is to satisfy these requirements, then since we 
must have, for any z, 


G=¢(«) = o[f(x)], 
it follows that ¢ must be so chosen that 


olf(x)J] = gl¢(x)] (7) 


for every x for which f(x) is defined. This condition is necessary and 
it is evidently also sufficient. We therefore inquire whether, for a 





es §¢ & ~~ 


= 
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given f(x) and a prescribed g(y), a function ¢(x) can be found sat- 
isfying (7). 

If the range of definition of g(y) is sufficiently great, such a 
function ¢ can always be found, and in infinitely many ways, as fol- 
lows. Let 2 be the absolute threshold intensity in the x-scale, and 
let yo be any positive number lying on the range of definition of g(y). 
Then define a sequence x; and a sequence y; by means of the following 
recursion formulas: 


Lie = f(%i), Yiu = 9(Yi). (8) 


Both sequences are properly monotonically increasing because of con- 
ditions (5). If at the upper limit of the x-range inequality (5) still 
holds, there will be a finite value, » — 1, of the index 7 such that 
f(%n-+1) equals or exceeds this limit. Then 2x, is defined to be this limit 
and the first of equations (8) is taken only for i < n —1. Then if 
the range of definition of g is great enough, or if y, is small enough, 
the sequence y; will include terms up to and including y, . 

If, on the other hand, at this upper limit for which f(x) is defined 
the inequality (5) must be replaced by an equality, then there will 
be infinitely many terms in the sequence 2; , and the limit of this se- 
quence is this upper limit, which we shall call x... Because the con- 
tinuity of f(#) we have f(x.) = x. In this case, if g satisfies (5) 
wherever it is defined, g must be defined for all positive values of y, 
and the terms of the sequence y; become infinite. But if g(y) = y 
for some value of y, then again the sequence y; has infinitely many 
terms which approach this value of y as a limit, whatever our choice 
of yo may have been. 

In either of the various cases, ¢(x) can be defined arbitrarily as 
a properly monotonically increasing continuous function of x over the 
interval 1, < x < 2, satisfying 


$(%) = Yo, $(4%,) =. (9) 


Now if x is in the interval x < x < «x,, and y is in the interval 
Yo <y < y:, then evidently x and 7 as defined by (3) and (4) satisfy 
%< 24 < %,and yi <9 < y.. Hence if ¢ is defined over the interval 
from x, to x. by means of the equation 


$(x) =9 (10) 


it is evident that (7) is satisfied for every x in the first interval 
from x to x,. We can now proceed sequentially to define 4 over the 
remaining intervals, using (3), (4) and (10) to define ¢ over any 
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interval x; < x < x;., after it has been defined over the preceding one, 
and equation (7) will be satisfied wherever it is meaningful. 

We have therefore shown that if there is given a single Weber 
function f(a) in terms of some scale units x, and if g(y) is an arbi- 
trary, preassigned function, both functions being properly monoton- 
ically increasing, continuous functions satisfying inequalities (5) ex- 
cept, perhaps, at the upper limits of the ranges over which they are 
defined, then it is always possible to determine, in infinitely many 
ways, a change of scale which transforms the Weber function f(x) 
into the Weber function g(y), provided only that the range over 
which g(y) is defined is sufficiently large. 

In particular, a single Weber function can always be given the 
Weber form by a suitable change of scale. 

It is reasonable to expect, however, that different individuals will 
have different Weber functions, though evidently for a given modality 
one would not wish to use different scales for different individuals. 
Thus presumably for one individual would be found a function f, in 
the x-scale, for another f, in the same x-scale, and if any transforma- 
tion ¢ is to be made to a y-scale, the same transformation would have 
to be made on both f, and f.. Given then, any finite set of empirical 
Weber functions, f, , f2, f;,---, will there always exist a scale trans- 
formation which transforms all these functions simultaneously into 
the Weber form (or some other specified forms 9, , 92, 93,°*:)? 

Before discussing this problem let us note a possible modification 
of the procedure for calculating the ¢ which transforms a single f 
into a single g. Instead of requiring that x, be the absolute threshold 
value, we could have taken as x, any value larger than this. For values 
of x greater than x, the procedure would be as before. In order to de- 
termine ¢ for values of x less than x, , we note that the equations (3) 
and (4) can always be solved uniquely for x and y as functions of 


x and ¥ respectively, because of the monotonicity and continuity of 
the functions of f and g. Hence along with the increasing sequences 
x; and y; defined by (8) we may place the decreasing sequences 2-; 
and y-; defined by 


ti=F(tia), Yyi=glyia). (11) 
Thus we have a sequence of intervals to the left of x. Now if 
G=9(2), to Sasa 


and if x and y are defined by (3) and (4), then 7, < x < 2, and 
¥..<y < y, and if we define y = ¢(x), condition (7) is again sat- 
isfied wherever it is meaningful. We may now proceed sequentially 
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over the intervals to the left until ¢ has been defined wherever f is 


defined. 
Hereafter we consider only functions g(y) of the form 


g(y) =ay, a = const., (12) 


and inquire only about transformations yielding the Weber form. Sup- 
pose we have two functions f,(x) and f.(x) , and suppose that for 
some x = x, on the range common to both functions we have 


f(x) > fev). (13) 


Let y be any positive number, and let the point (2, , yo) be our start- 
ing point for defining ¢. We have, now, two increasing «-sequences, 
x’; and x"; , defined by 


wi= fi (Xia), x"; = fe(Xi-). (14) 


Suppose that f,(x) is to be transformed into g.(y) = %y, and f(x) 
into g.(y) = ay. The two y-sequences are 

Yi = Yin, y", = BY" in, (15) 
and ¢ must be so defined that for every 7, 

yi = $(2#'i), y"; = o(2"5). 
Now from (13) it follows that x’; > x”; , and since ¢ is monotonic in- 
creasing we must have y'; > y";. But this is true if and only if 
a, > a. From this it is evident that the transformation cannot be 
made unless inequality (13) holds for all values of x on the common 
range. For if at some point x = 2) we should have the inequality in 
(18) replaced by equality, then we could use x, as our starting point 
and conclude that a, = a. 

If at some point of the range common to two empirical Weber 
functions f,(x) and f.(x) we have inequality (13) holding, then for f, 
and f, to be transformable to the Weber forms ay and ay , respec- 
tively, by the same scale transformation ¢(x), it is necessary that in- 
equality (13) hold over this entire common range, and that a, > a. 
This condition is by no means sufficient, as we can easily see. Thus, 
suppose the condition is satisfied, and, for simplicity, suppose 2» is 
the larger of the absolute thresholds. Then we have two sequences 
x; and x";. Because of (13) it is evident that not more than one 2’; 
falls at an endpoint of or within any interval (x"; , x”;..). Hence from 
the sequence x”; we can select a subsequence x”;, such that 


L"4, SUVS 2"; (16) 
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Corresponding to this subsequence is the subsequence y”;, = a,/*y, of 
y”; and evidently we must have 


Qnjiv <a," < ajvn, (17) 


since ¢ has to assume the values of the y’, when x has the value z’,, 
and also the value y”; when z has the value x”;. The above inequal- 
ities can be put into the simpler form 


_ loga, 


eo ia,” (18) 


jv<vacyprl, 





and this must hold for every positive integer » corresponding to which 
there is a term of the sequence x’,. The equality can occur in (18) 
only when the equality occurs in (16). If the equality occurs once, a 
is uniquely determined. If it occurs twice the two equations may or 
may not be consistent. And examples of sequences j, of integers can 
easily be found in which the inequalities (18) are not all consistent. 

If inequalities (18) are all consistent, then we may proceed as 
follows. Solve for x in all equations of the form 


x; = f(x), x; = f(z), (19) 


where x; is a member of one of the original sequences, or is one of 
the solutions resulting from the successive solutions of such equa- 
tions. If y; is the value of y associated with x;, then y;/a, or y;/d. 
is the value that must be associated with x. Proceeding thus we ob- 
tain additional points through which the graph of ¢(x”) must pass, 
and hence new inequalities in addition to (18) which a, and 4, ora, 
must satisfy. If all inequalities obtained in this way are consistent, 
and an a is chosen to satisfy them, then it may be possible to define 
¢ sequentially over intervals bounded by the points x and the points 
of the original sequences. 

When functions f, , f. ----- , f, are at hand, corresponding to ” 
individuals, one has initially sequences x; from the points of 
which one obtains other points by successive solutions of equations 
such as (19). From these one obtains certain inequalities on the pos- 
sible values of the a; aftey which one may be able to define ¢ as be- 
fore. Thus: 

If f,, f., ---- , f, aren empirical Weber functions for the same 
modality, corresponding to n individuals, it is not in general possible 
to determine a scale transformation ¢(x) in such a way that the trans- 
formed functions all have the Weber form. For this to be possible it 
is necessary but not sufficient, that, for every i and j, if f;(x) > 
f;(x) for some x, then this inequality holds for every x on the range 
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common to both functions. Further necessary conditions are obtained 
in the form of sets of inequalities by starting with a fixed x, , defin- 
ing additional points in all possible ways by x; = f(x.) where f is any 
one of the n functions, and either x; or x, is a point previously deter- 
mined. Corresponding to each x; so determined, the value y; that 
$(x;) must have can be determined in terms of an arbitrary initial 
y. and the constants a entering in the Weber forms. The inequalities 
arise from the required monotonicity of the function ¢. 

If f is supposed to have a particular analytical form as a func- 
tion of x and a single variable parameter, b , then for an individual 
whose parameter is b , the stimulus intensity x which is just notice- 
ably greater than zx is given by 


z=f(b,2). (20) 


In this case the question whether there exists a scale transformation 
y = ¢(x) such that f(b, x) takes the Weber form g = ay can be giv- 
en a very simple answer. For the equation 


y= wy (21) 


can be regarded as defining a transformation of the variable y into 
y, and we have in (21) a one-parameter continuous group of trans- 
formations. By a change of variable the group (21) is transformed 
into the group of translations 


z=ce+z2(z=logy, c=loga). (22) 


Now every change of variable transforms a one-parameter continu- 
ous group into a one-parameter continuous group, and conversely, 
every one-parameter continuous group can be transformed, by a suit- 
able change of variable, into the translation group (22) and hence 
into the “Weber” group (21) (see, e.g. 4, Vol. 2, Part 2, p. 89). Thus 
if f(b, x) is to be capable of taking the Weber form it must define a 
continuous group: 

The necessary and sufficient condition for a Weber function of 
the form (20), with continuous partial derivatives of the first order, 
to be transformable into the Weber form is that (20) defines a one- 
parameter continuous group of transformations. The scale trams- 
formation is, in this case, unique. 

Consider, now, a single individual, or a group of individuals re- 
garded as identical for purposes of the investigation, and suppose a 
large number of judgments are made comparing stimuli of different 
intensities. Suppose we do not fix a “discriminal” percentage, but in- 


quire, for each pair of stimuli, x and x, as to the probability distri- 
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bution of judgments “z > x”. For each x and for each x there will be 
a unique probability p of the judgment “x > x”, and conversely for 
each x and each p there will be an x such that the probability of the 
judgment “x > x” is. Thus we can express 


t=f(p,x), (23) 


where p is the probability of the judgment “x > x’, and x should be 
a properly monotonically increasing function of p for each x, and of 
x for each ». We shall suppose, further, that 


0=/f(0,z), o =f(1,2%) (24) 
x = f(4,2) (25) 

identically in x , and that, ruling out equality judgments, 
x=f(l1—p,2). (26) 


Since to every x and x there corresponds a p such that (23) is satis- 
fied, it is evident that, if we regard f as defining a transformation, 
the resultant of two transformations of the form (23) is again a 
transformation of the form (23). Moreover, these transformations 
include the identity (25) and the inverse (26). 

Hence, the function f(y, x) defines a one-parameter continuous 
group and is therefore reducible to the Weber form. 

Furthermore, since every one-parameter group is also reducible 
to the translation group, we can, by suitable change of variable and 
parameter, 


y = ¢(2) c=y(p) (27) 
reduce (23) to the form 
y=y-+e. (28) 


Now ¢ is a function of the probability of the judgment “x > x”, or 
“yj > y” , and c measures the “distance” 7 — y on the y-scale. This, 
however, is Thurstone’s psychophysical scale (11) in a somewhat 
more general form. Whether we think of introducing new scale units 
for measuring the physical intensity, or for measuring a psychological 
“discriminal process” as some functions of the physical intensity, the 
results are formally identical. 

Any probability distribution of judgments “x > x” when ea- 
pressed in the form (23) defines a one-parameter continuous group 
of transformations, and from such a distribution one can always de- 
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fine a psychophysical scale by (27) and (28), the scale units being 
expressed by an appropriate function y(p) of the probability p of the 


judgment “x > x”. 

So far as our acquaintance with the literature extends, the gen- 
eral question of the measure of the stimulus intensity is not ordinarily 
raised in discussing the “validity” of Weber’s law. Boring (1), how- 
ever, did raise this question, to which Johnson (7) replied that Fech- 
ner had clearly defined the stimulus as energy. Nevertheless, he ad- 
mitted that this cannot be, in general, and suggested that power 
should be taken instead. But whether it be power or energy or any 
other, is this to be measured at the periphery, at the first-order neu- 
rons, or elsewhere? 

Cobb (2) points out that the form of the Weber function varies 
with the setting in which the stimulus is placed and remarks that 
while for one setting the Weber law may be verified over a wide range, 
it will not necessarily be verified when the setting is changed. But 
from the interpretation of the “Weber law” given here, this conclu- 
sion is not obvious. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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SCORING WEIGHTED MULTIPLE KEYED TESTS 
ON THE IBM COUNTING SORTER* 


RALPH BEDELL 
THE UNIVERSITY OF NEBRASKA 


Tests or personal inventories with differential item response 
weights may be scored by means of punch card equipment. Detailed 
instructions are given for preparing the cards and scoring the forms. 
The scoring speed is approximately four to eight times that attained 
by manual scoring. 


The method of scoring herein described for weighted multiple 
keyed tests by use of the International Business Machines Counting 
Sorter can be used by any clerk trained to do sorting and counting on 
the machine. The machine has a relatively low rental cost and is 
usually used for a wide variety of statistical work other than test 
scoring. This type of test scoring does not require modification of the 
machine nor in the administration of the tests. While the method 
may be adapted to any multiple keyed test such as an interest inven- 
tory or personality questionnaire, detailed explanation is here given 
for scoring the Vocational Interest Blank for Women by E. K. Strang. 
This test can be completely scored for all eighteen keys by one clerk 
at the rate of two tests per hour by using a machine operating at 220 
cards per minute. Two clerks using one machine operating at 440 
cards per minute can score for all eighteen keys and profile the re- 
sults at the rate of four tests per hour. This scoring spee dis approxi- 
mately that of other machine methods now being used, and four to 
eight times that attained by manual scoring. The scoring speed is 
not affected by the number of tests to be scored; consequently it is 
as useful in clinics scoring small numbers of tests at a time as in mass 
testing programs. 

The method consists of punching one Hollerith card for each test 
item. The punches are coded to give the number of the item, and the 
weight for each possible response to that item in each key for the test. 
The cards are then sorted by hand according to the testee’s responses 
to the test, and the various weights totaled for each key on the IBM 
counter. Totals are entered on a:tabulation sheet on which the re- 
quired simple computations are made. 


* The author is indebted to H. M. Cox, University of Nebraska, for certain 
suggestions on the columnar coding and the method of checking herein presented. 
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Coding 

To code the weights for each response to each item, it is necessary 
to provide for nine weights varying from — 4 through +4. Only four 
positions on a column are required to code the nine weights if com- 
binations of +1, —1, +3, and —3 are used. These weight code equiv- 
alents with the Strong weight given first are: +1 = +1;-+2 = —1, 
+3;-+3 = +3;+4 =+1,+3;0 = no punch; —1 = —1;—2 = +1, 
—3 ;—3 =—3; —4=— 1, —3. Since the card has twelve possible 
positions in each column, three keys can be coded in one column for 
one response. The three possible responses require three columns for 
these keys. Each succeeding three keys are coded in each succeeding 
three columns. Thus, for the eighteen keys, eighteen columns are 
needed, plus three columns for the number of the item. 

To make it possible to construct a convenient record and tabula- 
tion sheet to be used by the operator, it is desirable to code the three 
keys in each column as follows: 


Position Weight 
on Card Code 
1 +1 
(Key A) 2 1 
3 +3 
4 —3 
5 we | 
(Key B) 6 —1 
7 +3 
8 —3 
9 +1 
(Key C) 10 (0) —1 
11 (x) +3 
12 (y) —s 


The accompanying tabulation sheet (Figure II) shows the meth- 
od of coding for the eighteen keys for this test. The first three col- 
umns are reserved for the number of the item. Column 4 carries the 
weights for the L response in the keys for Artist, Stenographer-Sec- 
retary, and English Teacher; column 5 carries the weights for the I 
responses in these keys; and column 6 the weights for the D responses. 
The system is illustrated by the coding for item 109 Stenographer 
L. I. D. The Strong weights in each key for each response with the 
columns in which these weights are coded are given in Table I. The 
complete punching of the card for this item is given with the column 
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numbers in bold-face type and the position numbers in light face. 

1—1; 20; 3—-9; 42, 4, 5, 7, 0; 5—2, 6; 6—2, 3, 6, 8; 7—6, 0; 

8—O; 9—2, 5, 9; 10—2, 6; 11—1, 5; 12—no punch; 13—1, 4, 5, 8, x; 

14—-2; 15—2, 3, 5, y; 16—5, 0; 17—5, 0; 18—5, 8, 0, x; 19—2, 6, 
0; 20—1, 5; 21—9. This punched card is shown in Figure I. 


/ 109 ‘ 


i i 
DP 
| PRRRRREEE! BF DEERE) DROP RE ER RREO ROR RERRReREEEE 
2220022002 228222222222 22222222222222222222222 
33333933333333933333333333333333333333333333333333 
4440444444469 44446444644644444464464446444464646446444644444 
SSSMSSSSMSMSMSHRURSH555555555555555555555555555555 
GECCMBRGCCHGGCC6C CEC MPECEECEC6EGE6GCEG6GEGE6GSEG6GGEEGGEE 
7797777797777999797979779979799777797777979777777977997777 
BEEK PSSESCSSHSSSSPSSESCSSSSSESEKKEBSEKBSRKESEEEEBEE 













999999999999999999999999999999999999999999999999 
2 Fess hae" THANE YONPELONAAHASHUN NADH NSORnennesees 








FIGURE I 


The completely punched card for item 109. The position next above 0 is x and 
the top position is y. 

One card is punched for each item making a deck of 410 cards. 
This deck is then interpreted* on the first three columns so that the 
number of each item appears on the card above columns 1, 2, and 38. 
The cards are then duplicated* and verified by the Reproducing 
Punch* and interpreted* to make as many decks as are needed. To 
avoid possible loss and excessive wear the deck of master cards is 
kept in a safe, and removed only to make the duplicate sets which 


are used for scoring. 


* These are card-making operations, not test-scoring operations. The Inter- 
preter is a machine which automatically prints at the top of each card informa- 
tion punched in the card. The Reproducing Punch is another machine which au- 
tomatically duplicates and verifies new decks from the original punched deck. It 
is not necessary to lease these machines as the service here indicated may be ob- 
tained through any International Business Machines office at nominal charge. 
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TABLE I 


Original Strong Weights in Each Key for Each Response to Item 109 with 
the Columns in Which the Weights Are Coded 
















































































Column | 4 5 6|— . . - {10 1 RB 
Response | L I D . * = LT = 
oe ee one ail ~ | Teacher 

Artist |—4 —1 +2 | Lawyer 0 0 —1/inGeneral |—1 +1 0 

Stenog’er- | > oe | Math.-Phys. 

Secretary +4 —1 —4 ] Dentist —1 0 +41 | Sei. Tchr. —1 +1 0 

English —<—- > ee | Insurance | Femininity- Fel 

Teacher |—1 0 0 | Saleswmn. |—1 —1 +41 | Masculinity 0 0 0 

Column | 13 14 15 16 17 18 19 20° 21 

Response | L I D| a as L IB 

» Sat te ee es aaa Soc. Science 

Author  (|—2 —1 +2/]Nurse | 0 0 0 Teacher —1 +1 0 
Y.W.C.A. | 

Librarian | —2 0 +1 | Housewife | +1 +1 —2 | Secretary ssl +1 0 

Ofice | =~ ~~—~—~S~*~CSocial’ == " 

Worker +3 0 —3 | Worker —1 —1 +2) Physician —1 0 +41 








Table reads: Code in column 4 the L responses for Artist, Stenographer-Secre- 
tary, and English Teacher; the original Strong weights for which are —4, +4, 
and —1__sirespectively. 


Checking the Card Punching 

While the accuracy of the original card punching may be checked 
in the usual manner on the verifier punch, a trial balance technique 
for using the counter to check accuracy has been developed. A bal- 
ance sheet is set up which carries the algebraic totals of the weights 
provided by Strong for each twelve items for each response for each 
key. This gives thirty-six totals for each twelve items. These totals 
are then cumulated for each sixty items. The cards are then run 
through the counter in groups of sixty for each response for each key. 
If totals do not check, an error is indicated and that group of sixty 
is rerun through the machine in groups of twelve to check for these 
totals. The group of twelve which contains the error is then read vis- 
ually for the key and the response on which the error is indicated. 


Scoring Procedures 
The scoring procedures are given in detail for one clerk using a 
machine running at 220 cards per minute. 
1. Record on the tabulation sheet (Figure II) the name of the 
person whose test is being scored. 
2. Hand sort cards into three stacks corresponding to the tes- 
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FIGURE II 
Tabulation Sheet 
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Name 
Directions: Record total for each counter opposite its position number and under the columns 


on which the count was le. 


8; equals the algebraic sum of the numbers opposite 3 and 4, or 7 and 8, or 11 and 12. 


381 equals three times 41. 


8 equals the algebraic sum of the numbers opposite 1 and 2, or 5 and 6, or 9 and 10. 
Tot. equals the algebraic sum of 35: and 8&2. 


Record 


the correct signs for 8:1, 381, S2, and Tot. 
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This Tabulation Sheet is expanded and mimeographed to fill an 814 x 11 sheet. 


tee’s responses. This places all L responses in the first stack, 
all J responses in the second stack, and all D responses in the 


third stack. 


3. Count the L stack on column 4 and remove this stack from 


Op 


the machine, keeping it separate from the other two stacks. 
Without resetting the counters, count the J stack on column 
5 and remove this stack as above. Count the D stack on col- 
umn 6 in the same manner. 

Record the counter totals on the tabulation sheet. 


Reset counters to zero. 
Repeat 3, 4, and 5, except count the L stack on column 7, 
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the J stack on column 8, and the D stack on column 9. Con- 
tinue this procedure, making sure that totals are recorded 
and counters reset to zero after the counts are completed 
on columns 6, 9, 12, 15, 18, and 21. 

7. While the machine is making the above counts, perform the 

computations indicated on the tabulation sheet. 
8. After all counts are completed, place all cards in the machine 
and resort to the original numerical order by sorting on col- 
umns 3, 2, and 1. 

9. While performing operation 8, note and hand sort the sec- 
ond deck as indicated in 1 and 2. 

For a machine operatinz at 440 cards per minute, two clerks are 
necessary to keep the machine running with minimum time loss. One 
clerk operates the machine and records the totals. The other clerk 
does the hand-sorting and computations indicated in 1, 2, and 7 
above, and in addition profiles the results or writes them on the test 
blanks according to the system desired. 


Some Precautions 

To save time due to wear on the cards, it is desirable to have 
available five decks for each machine. One of these decks is set aside 
to be “robbed” of individual cards that may be unduly roughed or 
torn in the machine. Replacing such worn cards in the deck being 
used is much more economical than trying to force worn cards through 
the machine or stopping to reproduce the cards on the Alphabetic 
Duplicating Punch. Such removed worn cards can be set aside, later 
duplicated, and the new cards returned to the “robbed” deck. As 
clerks sometimes find time to keep their sorting ahead of the machine, 
the remaining four decks are used for scoring. Clerks should handle 
all cards by the ends as handling on the sides increases machine wear. 
New cards should be examined to make sure they have the proper 
bend for the machine before being put into use. No errors made by 
the IBM Reproducing Punch have been found, and it is believed that 
checks of card punching beyond those indicated for the original deck 
are unnecessary. The accuracy of the counters on the machine can be 
checked by counting all of the cards on any one column, the totals for 
which have been determined previously when the machine was known 


to be accurate. 


Adaptability to Other Multiple Keyed Tests 
While the method has been given in detail for only the Strong 
Vocational Interest Blank for Women, adaptation can readily be made 
to several other tests. The Strong Vocational Interest Blank for Men, 
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the Humm-Wadsworth Temperament Scale, and the Bernreuter Per- 
sonality Inventory illustrate the types of tests for which this scoring 
method is appropriate. For any test which has the same range of 
weights as the Strong blanks, it is believed that the above coding is 
the most economical that can be derived for this method of scoring. 
Tests that have a different range of weights will require different cod- 
ings, but the fundamental principle remains the same. The scoring 
speeds mentioned at the beginning of this article can be compared 
with the speeds of other methods being used to determine in advance 
whether or not it would be economical to use the method here pro- 
posed. By checking the time clerks require for the hand-sorting, com- 
puting the time required by the machine, and estimating the extent 
to which operations may be carried on simultaneously, it is possible 
to arrive at a reasonably accurate estimate of the scoring speed that 
may be attained for any test. 
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CONTRIBUTIONS TO THE MATHEMATICAL THEORY OF 
HUMAN RELATIONS Iil 
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In continuation of a previous study, somewhat more complex re- 
lations are considered for the case of interaction of two classes, of 
which one “organizes” the activities of the other. By way of illus- 
tration it is shown how different assumptions concerning the atti- 
tudes of the individuals to such an interaction of classes lead to 
different mathematical expressions. 


In several previous publications (1, 2, 3, 4) we have attempted 
to describe in mathematical terms some social phenomena. As has 
been emphasized in those papers, such a description is at present pos- 
sible only at the expense of a very great oversimplification and ideal- 
ization of the phenomena to be described. As in every theoretical sci- 
ence in its early stages, we must confine ourselves to rather simple 
conceptual situations, which may bear only a very remote resem- 
blance to actually observed phenomena. We hope eventually to come 
to the study of sufficiently complex situations, so that our formulae 
may be applied to real cases. But it is impossible to arrive at this 
result without a previous purely abstract study, that is without first 
studying the theory for theory’s own sake. A more thorough discus- 
sion of this assertion is to be found in the above mentioned previous 
publications, as well as in other papers, to which references are given 
in the above mentioned ones. 

The present paper is a continuation of the previous ones, and 
therefore presupposes a knowledge of these. In particular we shall 
have to refer to results obtained in the last two papers (3, 4), which, 
for abbreviation, will be referred to as HR I and HR II respectively. 

In HR II we have investigated the interreaction of two types of 
individuals, who are characterized by the following purely hypotheti- 
cal properties: 

Each individual of type I produces per unit of time an amount 
p, of goods necessary for the maintenance of his life, and he consumes 
per unit time an amount ¢, of these goods. Similarly every individual 
of type II produces and consumes correspondingly amounts p, and 
¢,. Denoting by N, and N, the numbers of individuals of type I and 
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II respectively let us consider the case, in which 
Ni.<< N23 Mm—A =a >0; P—e=e<O0. (1) 


Then in the absence of any exchange or external supply, the individ- 
uals of type I are capable of indefinite existence, those of type II 
are not. 

Now let the individuals of type I be capable of organizing their 
work, as well as the work of individuals of type II in such a way that 
the production of necessary goods is largely increased. In other words 
let the production of goods require a certain amount of organization 
skill, which is possessed by individuals of type I but not by those of 
type II, so that when left to themselves the individuals of type II 
cannot produce enough goods, but they can do so under the direction 
of the individuals of type I. 

In HR II we have investigated the following relation between 
the individuals of type I and II: The latter agree to work under the 
direction of the former, so as to increase their output of goods, pro- 
vided they do receive a part of the goods thus produced for their con- 
sumption. On the other hand, the individuals of type I will agree to 
direct the work of the individuals of type II only if they receive also 
part of the goods thus produced. For simplicity’s sake we have con- 
sidered the situation, in which every individual of type II, who ex- 
pends an amount of labor w. per unit time, gets a fraction éw. of 
goods per unit time. The quantity 6 may be considered as the price of 
his labor, paid not in money, but in goods produced by himself. 

We shall consider in this paper a somewhat more complex and 
perhaps more natural situation, in which an individual of type II re- 
ceives goods not only in proportion to the amount of labor he expends 
per unit time, but also in proportion to the amount of goods which he 
produces per unit time under the direction of individuals of type I. 

In HR II we considered that the total amount of goods produced 
per unit time by all individuals of type II under the direction of the 
individuals of type I is equal to (HR II, equ. 2): 


P=AN.Ww,, (2) 


where A is a function of the ration » = N,/N., A = f(y), and f(m) 
has the following properties: 


f(m) >Ofory>0; f(0) =0, fw) =const. (3) 


The reasons for those properties have been explained in HR II. 
Expression (2) tacitly implies, that the amount of supervising 
work done by individuals of type I per unit time is constant. It is of 
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interest to consider the more general case, in which P is proportional 
to the product w,w, of the amounts of labor w, and w, given per unit 
time by individuals of type I and II. P by definition is the amount of 


. goods produced per unit time by individuals of type II under the di- 


rection of individuals of type I, or in other words P is the excess over 
the “natural production” N.p, of the individuals of type II when left 
to themselves. Therefore P should vanish for either w, = 0 or w. = 0. 
Thus we shall now consider the case, where 

P= f( 7" Nz WiWe , (4) 


2 





f(y) being the same as before. 

By an argument similar to that used in HR II, we now find, in- 
stead of equations (10) and (11) of HR II, the following equations 
for the rates of change of the total amount of goods W, and W, pos- 
sessed by all individuals of type I and II correspondingly: 


aw, 


ae Nes + No F(a) WW (1 — 6) ’ (5) 
baad = Nye. + ON2 f(y) wie . (6) 


From the definition of 6 we always have @ < 1. 
In HR II we have supposed that the “price” 6 is related to the 
amount of labor w, by a simple hypothetical demand equation: 


be 
Wz = Woe wees (7) 


where w,, and b, are constants. A similar situation may be considered 
for w,. Individuals of type II will be willing to give more supervis- 
ing labor, the greater the amount of additional goods they receive in 
return. We shall however introduce here in our assumptions an ex- 
plicit asymmetry in the behaviour of the two classes. This is partly 
made in order to simplify some calculations. The effects of not intro- 
ducing such an asymmetry will be discussed at the end of this paper 

We shall assume that while the individuals of class II determine 
the amount of labor they are willing to give by the fraction of the 
goods which they receive in return per unit of goods produced, the 
individuals of class I determine the amount of supervising labor which 
they are willing to give by the total amount of goods, which they re- 
ceive from class II per individual of class I as the result of coopera- 
tion. Thus we shall put 
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b; 


un 9 PU ON; * 





Introducing for P its expression (4), we find: 


bin 


7(q)w,w.(1— 8) * (8) 


W; = Wo — 





Equation (8) may be written: 
f(y) (1 — 0) ww? — Wo (y) A — 6) ww, + O47 = 0. (9) 
Solved with respect to W, , it gives: 





Wor + V Ww 01 f?(m) (1 — 6)? w*, — Abi nf () (1 — @) wm, : (10) 
y) 2 f(y) (1— 4) w, 


We shall now discuss the second term on the right side of equa- 
tion (5). This term represents the amount of goods received per unit 


time by all individuals of type I from those of type II. To this end 
consider the expression w,w.(1— 6). From equation (10) and (7) we 


find after rearrangements: 





wW,= 














WW. (1 — 6) —_ at ae 6) (Woe — oy 
| rat Abin 
Mit ese ey ee 
The expression 
x= (1— 86) (Wo. — 02/8) , (12) 


which occurs in the denominator under the radical sign has a maxi- 
mum for 








0a be 
. (13) 
The maximum value 2.2 of « is equal to: 
Lmaz — Woe2 — b, me 2 V b2Wo2 . (14) 
If 
4bin 
Tn) Watmar <1 
or 
4b.n 


ma > SI 15) 
" f(y) wo ( 
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then the expression w,w.(1— 96) has two real positive values for a 
given @, provided @ lies within an interval 


A<0<a<1. (16) 


which is in the neighborhood of the value given by equation (13); 4, 
and 6, are roots of the equation obtained by equating the right hand 
side of equation (12) to the right hand side of inequality (15). 

If inequality (15) does not hold, then w,w,(1— 9), and there- 
fore the last term of equation (5) is imaginary. Under these condi- 
tions a cooperation of the two groups is impossible: no agreement as 
to price can be reached. It is interesting to note, that inequality (15), 
which represents the conditions for the possibility of cooperation, in- 
volve such quantities as the function f (7), or the constant w»: , which 
are more or less of a physical nature, as well as the constant b, , which 
is more of a psychological nature. 

The existence of two values of the expression w,w.(1 — 9) is 
due to a two-valued nature of w,, as defined by equation (8). This 


= bn 
F(n)w, Wo C1 - @) 








WwW, 


FIGURE 1 


however is not connected with the specific analytical expression chos- 
en for w,, but lies rather in the nature of the general assumption 
made, namely that the relation of w, to the return in goods is such, 
that for returns below a certain value the individual will not wish to 
give any labor, while for increasing return the amount of labor in- 
vested increases, tending to a limiting value determined by the phys- 
ically possible maximum of the labor given. This, combined with the 
assumption that the return in goods is in its turn proportional to the 
amount of labor given, leads to the two valued character of w, in 
terms of 6. For the left hand side of equation (8) is represented by 
a straight line (Figure 1), while the right hand side is represented 
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by a curve, which, regardless of its detailed analytic specification, in 
general, intersects the straight line in two points. 

Returning now to equation (11), we see that the value of 
W,W.(1 — 6) which corresponds to the upper sign before the radical 
is always the largest of the two. For 6 = 6, as given by equation (13), 
both the expression in brackets and the expression in front of it have 
a maximum. Hence for 6 = 6,, , the quantity w,w.(1— 6) has a max- 
imum for the upper sign before the radical. For the lower sign it 
may or may not have a maximum. 

In HR I we have established some equations describing the con- 
ditions under which one of the two classes controls the behaviour of 
the other regardless of such interreactions as studied here. If the re- 
lations between the constants, involving social control (HR I), are 
such, that in the present case the class I, which is composed of indi- 
viduals of type I, controls behaviour of class II, (which is not a priori 
necessary), then class I will prescribe as price for the labor of class 
II the value 6,,, for this makes dW,/dt a maximum. 

Due to inequalities (1), dW,/dt is always positive for 6< 1. But 
dW./dt may be negative. If 6 = 6,, then, in order to make dW,/dt 
also positive, we must have 


N2eé2 
Nf (m) ’ 


Otherwise the “wealth” of class I will increase while that of class II 
will decrease. Class II will not be able to exist indefinitely, and with- 
out it class I will not be able to have additional goods also. For a con- 
tinued negativity of dW./dt will result in a decrease of N., and this 
will result in a decrease of P and of the last term of equation (5). 
Thus condition (16) is necessary, but not sufficient for the coexis- 
tence of such two classes. It must be supplemented by condition (17). 

If class I fixes 6 at 6,,, this will result in a maximum of dW,/dt. 
But this does not mean that the rate of increase of the total goods for 
both classes, that is d(W, + W.) /dt will be a maximum. In fact, this 
will not be the case. From equations (5) and (6) it follows that 
d(W, + W.) /dt has a maximum for such a value 6’, of 6, for which 
W,W, has a maximum. Denote 


Y.=WW.2>03 Y= ww.(1_—s) =y,(1—0) >0. (18) 


The quantity 6’, is defined as that value of 6 which makes y, a maxi- 
mum. Hence 


6w,w,. > — ford = ¢,, . (17) 


dy; 


>, = 6 for 0= 6'n . (19) 
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dy, _ ‘i ey 
If 6,, is the value, which makes y, = y,(1— 6) a maximum, then 
dy, _ = 
“aa =0 for 6 “ae 
But for 6 = 6’, , because of (19) 
dy. a dy, _ 
Fe 91 (8m) - CL — Om) a6 =—y,<0. (20) 
Hence @'m > Om. 


In other words, under such conditions as are described by equa- 
tions (5), (6), (7) and (8), such a cooperation of the two classes, 
which is optimal for the population as a whole will not be optimal for 
class I. 

Denoting by @”» the value of 6, which makes 6w,w, a maximum, 
and therefore gives a maximum dW,/dt, we can prove by a similar 
argument that 


wa > Fs > Am ° (21) 


If the class II controls the behaviour of the population, or if a 
third class is present, which, while not participating directly in the 
cooperation discussed here, influences nevertheless the behaviour of 
the population by virtue of a social mechanism discussed in HR I, 
then 6 may be fixed either at 6’, or at 6”,,. This will lead to optimal 
conditions for either the whole population or for class II. It may how- 
ever happen, that 6”, > 6, or even @’,, > 6., [Equation (16) ]. In this 
case no cooperation will be possible at all. Class I will continue to 
exist due to «, > 0, though the rate of increase of its wealth will be 
very much reduced. Class II will not be able to exist indefinitely. Un- 
der the conditions studied here, the coexistence of the two classes may 
be possible, for a certain choice of constants, only when it corre- 
sponds to conditions close to those optimal for class I. 

If we do not introduce the asymmetry, expressed in the different 
structure of equations (7) and (8), but take for w, an expression 
similar to (7), except that 6 would be replaced by (1 — 6), we would 
have simplified the situation somewhat. We now would have: 


(1— 6) w,t, = (1— 8) (20 — >">) (wor — "2, (22) 


1—9@ 


which has a maximum for 





Wo1Wo2 


6=06, = .| —— + 
(Wo, — 5,) b. 


(23) 
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Similarly as before we shall find 6”’» > 6m > Om. 

On the other hand, if we take for both w, and w. expressions of 
the form of equation (8), then the situation is mathematically more 
complex. Though essentially similar results are obtained, yet closed 
expressions are very cumbersome, for we have to deal now with two 
simultaneous quadratic equations, which reduce to a single fourth de- 
gree equation. 

The different expressions that may be used for w, and w, corre- 
spond to different assumptions about the psychological attitudes taken 
by each class towards cooperation of the type described. The relations 
discussed here are too simple to be of any practical use. They merely 
illustrate how different psychological attitudes of different social 
classes may be in principle translated into mathematical language. 
But it is only through the study of such simple illustrations that we 
may hope to arrive at the study of more concrete real cases. 

The author is indebted to Dr. A. S. Householder for a discussion 
of the paper and valuable critical comments. 
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CORRELATION COEFFICIENTS FROM THE 
FACTORIAL MATRIX* 


P. S. DWYER 
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This paper shows how to compute multiple correlation coeffi- 
cients, partial correlation coefficients, and regression coefficients 
from the factorial matrix. Special emphasis is given to computa- 
tion technique and to approximation formulas. The method is ex- 
tremely flexible in application since it may be applied to any subset 
of the original set of observed variables. It is also extremely useful 
when many of these coefficients are desired. 


PART I 


MULTIPLE CORRELATION COEFFICIENTS FROM THE FACTORIAL MATRIX 

1. Introduction. One of the disadvantages of multiple correla- 
tion theory is the enormous amount of time needed to compute the 
multiple correlations when nv is large. Even the Doolittle method, 
(1, 2, 3) which is being recognized as an efficient method of comput- 
ing multiple correlations, demands a large amount of computation 
when v is large. It is shown in this part how multiple correlations 
can be obtained from the factorial matrix, and in particular, from a 
centroid solution, though any solution which gives communality and 
uniqueness factors is satisfactory. Though a centroid solution does 
take time, the total time necessary to obtain the multiple correlations 
is but a small fraction of that demanded by other methods if the num- 
ber of variables is large, if the number of factors is small, and if a 
number of such coefficients is desired. 

2. Notation. (a) From multiple correlation theory. We assume 
that there are n standard variables, x, , %2,-++ , Vi,+++, Uj, ++, Un. 


*Editor’s Note: The reader will recognize the subject-matter treated in this 
article as closely allied to that of an article by Louis Guttman, “Multiple Recti- 
linear Prediction and the Resolution into Components,” in the June, 1940 issue 
of this journal. Although there is considerable overlapping in the topics consid- 
ered in these two articles, it was felt that some readers would be more interested 
in one approach, while others would gain more from the other approach and that 
still others would find both presentations of value. Guttman’s article is somewhat 
more concerned with the theoretical aspects of the problem, while Dwyer’s article 
emphasizes the technique of computation and approximation. The manuscript 
for Dwyer’s article was received while Guttman’s article was in press. At Dwyer’s 
request, Guttman made available to him a pre-publication copy of his manuscript 
so that Dwyer was able to insert in his manuscript references to the earlier one. 
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The variables are recorded for a population of N individuals and the 
matrix of the correlations, [7;;], is computed The determinant of 
this matrix is indicated by A and the determinant of the matrix 
formed by deleting variable x; is denoted by 4;. If R; indicates the 


multiple correlation of x; with the n — 1 remaining variables, then 
by the usual formula (4, p. 301) 
vidal 
Rg, = if 1— a (1) 


Ri, is used to indicate the multiple correlation of x; with all the re- 
maining variables except x;. Similarly R;)i,, indicates the multiple 
correlation of x; with all the remaining variables except 2; and 2;. 

(b). From multiple factor theory. We assume that the corre- 
lation matrix has been subjected to an analysis which results in 
r + n standard orthogonal variables, r of which are common vari- 
ables and » of which are uniqueness variables. More specifically we 
assume that every x; can be represented by 





Lj = AjyYs + VjoYo + +++ AjeYR + +++ + AY, + UjZ;, (2) 
where ¥; , Yo,°**, Yr, 215 22, °**, Zn are the orthogonal variables, and 
the values a; are known. 

It follows at once, since 7;; = oe , , that 


ij = UjyQj, + AjgAjg + +++ AiO, +--+ 4;,0;,, tI (3) 
and 
L = jy? 4 Wja? +++ Ayn? + ++ + Gj? + U;?, (4) 
as indicated by Thurstone (5, p. 66). Furthermore we define h;? = 
@;,7 + ajo? + ----+ a;,? to be the communality of variable 7 and u;* to 
be the “uniqueness.” Then 
h?+uz=1 and wu = V1i—h;. (5) 


We next define the common part of «x; to be that part of x; which 
is expressible in terms of the common variables y,. Then 


5 = AY + Aj2Y2 + +++ DixYu + +++ VY. 


It is to be noted that x’; is not a standard variable since its variance 
is not unity. However, if every factor loading a; is divided by h; we 
get the standard variable 


Tr 
2", = = A jx ’ 
=1 
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where Aj, = ~ . The factor loadings are increased by this process 
J 

and hence we use the term “augmented” (5, p. 158). Furthermore, 

we introduce the term “communality augmented” to distinguish these 

from the “uniqueness augmented” loadings introduced in this paper. 

We let U? = [u;?] indicate the diagonal matrix composed of the 

uniqueness loadings. Then (3), (4), and (5) can be combined to give 


the fundamental matrix equation 
R =AA'+ U?. (6) 


8. The value of R;. R; can be expressed in terms of the factor 
loadings if A and A; can be so expressed. Now 4A; can be so expressed 
if A can. Hence such an expression for A will, essentially, solve our 
problem. A general development using matrix terminology is pre- 
ceded by a special case, n = 3 and r = 2, in which determinant the- 
ory is used. We first express 4 in the form 


| 1 Ti2 113° yi =Ane 
| 


| Tr L Tes Ger Aes 
A= | Ts1 Tso 1 Qs, Ase (7) 
| . 6¢ 8 2 86 
| 4; ea. 2 | 


and eliminate the correlation coefficients by multiplying column 4 by 
—t,,, column 5 by —a,., and adding to column 1; by multiplying 
column 4 by —a.,, column 5 by —a.., and adding to column 2; by 
multiplying column 4 by — a;, , column 5 by —a;. , and adding to col- 
umn 3. We then have 





u,? 0 0 G@ Ge | 
Us? 0 = Aa, Age | 
A= | 0 0 ‘* Gy Ge | ‘ (8) 
—O,, —M, —d;, 1 0 | 
—t4 —te te 0 1 | 
We now divide the first row and the first column by u,, the second 


row and the second column by w,, the third row and the third column 


a; 
by u, , and we get, if we let — = B;, , the “uniqueness augmented” 
j 


loading, 





¢ 
% 
f, 
% 
' 
’ 
BA 
4 
a 
4 
Hy 
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1 0 A te BK 
0 1 = s * 
0 0 PP oe 

eg ally lly Fd 

eee OE ee ee | 


This can be further simplified by multiplying the first row by 
B,, , the second by B., , the third by B;, and adding to the fourth row; 
by multiplying the first row by B,., the second by B.,, the third by 
B;. and adding to the fifth row, so that we have 


1 + SB;,’ SB iiBic 
A = U,2Uy?U3? : (10) 
SBiBi2 1+ SB? 


A = Uy? U2? Ug” 


(9) 





This can be written as 


M,, Ma 
A = u4,7U,*u,? Ms Me . (11) 
In the general case we have 
A{R] = 4 lor |=4 meme & =A [ar | - aqwi4| pr 
=A 14] . |= 4 [w2] A [+ BB’). (12) 
+ BB 


The determinant, A[J + BB’], is a determinant of 7 rows and r col- 
umns, the elements of which are 


| 3 BuBu if kl | 
Mu=)"™ , Le (13) 
| 1+ > BiBi if =f 


The matrix [J + BB’] is the equivalent of Guttman’s matrix Q (12, 
section 4). We may use the above result to indicate the value of 
A; if variable 7 is deleted. Then (1) becomes 


a , Atl + BB 
R= || Lu) a (14) 





where 4;(J + BB’) is a determinant having r rows and r columns 
whose elements are given by 
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S BiB —BrBn if kl 
Mi) 5 = Se ws ° (15) 
1+ > BuBi — BrBi if ee 


See also (12, Section 4). 

The formula (14) is easy to apply (particularly if r is small) as 
each element of A; can be obtained from the corresponding element 
of A by a simple subtraction of a product. The flexibility of the for- 
mula is also apparent since the correlation of any variable with any 
subset of the original variables is easily effected. It is only neces- 
sary to complete the values of A and A; from the factor loadings of 
the variables included. 

4. Illustration. As a first example we use an illustration of 
Thurstone (5, p. 128) in which the six-place residuals are 0. In the 
factorial solution the correlation matrix involving 7 variables is re- 
placed by a factorial matrix of two factors. The factor loadings, com- 
munality, and uniqueness for each variable are given in columns (2), 
(3), (4), and (5), while the “uniqueness augmented” loadings are 
given in columns (7) and (8). 


TABLE I 
Computation of “Uniqueness Augmented” Factor Loadings 
(6) (7) 


(1) (2) (3) (4) (5) (8) 


Variable 

j a; dio h;? u;? u; B;, Bj, 

1 .659828 -120945 .450000 .550000 .741620 .889712 -163082 
2 .830332 -265611 .760000 .240000 .489898 1.694908 542176 
38 —.541290 -637969 .700000 .300000 .547723 —.988255 1.164766 
4 —.126124 -770774 .610000 .390000 .624500 —.201960 1.234226 
5 437356 -590526 .540000 .460000 .678233 644846 870683 
6 .687638 —.336776 .520000 .480000 .692820 .920351 —-.486094 
7 .904489 109084 .830000 .170000 .412311 2.193706 .264567 


In getting the value R,, we find that 
M,, = 11.756955, Mi. = M., = .858146, M.. = 5.264916, 
Mj;)1, = 10.965368, Mys)1¢ = Moi). = .218050, 
Mos). = 5.238320, 
so that A[1 + BB'] = 61.771112, A, [J + BB’] = 57.394716, 


A[I + BB’] 


and 





R, = V1— (.550000) (1.076251) = .638797. 
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This multiple correlation coefficient would be tedious to compute 
by conventional methods, so, in order to provide a check which the 
reader can follow with little computation, I select the case in which 
x, is predicted from x, and x;. In this case 


M,, = 5.640949 , My. = M2, = —.087051 , M.. = 2.677230, 
Mas, = 4.849861, Myo), = Mory.¢ = —.232147, Moo). = 2.650635, 
with A[I + BB’] = 15.094540, 4,[{J + BB’] = 12.799994, and 





R,(= Ry23) = V1 — (.550000) (1.179261) = .592795. 





Ay. A; 
—=F's3 


» we get 
an : “ih 








If we use the formula R, = J Tes 





bey (.58) (.5016) + (.28) (.1176) o 
R, 5216 = 592795. 
The conventional formulas are preferable for the case of three vari- 
ables, but the advantage of the new method increases as the number 
of variables increases. For example, it was possible to obtain the 
seven multiple correlations obtained by estimating each variable from 
all the others from the three-place factor loadings in about three- 
quarters of an hour. The work is outlined in Table II. 

There are n(2"* — n) multiple correlation coefficients. Any one 
of these can be obtained from the uniqueness augmented factor load- 
ings. Thus in the illustration there are 399 different multiple corre- 
lations, any one of which can be obtained in a relatively short time. 

5. Approximate factorial solutions. The problem above is one 
in which the residuals are 0, to the required degree of approxima- 
tion, and the multiple correlations agree to the desired number of 
places with the multiple correlations obtained by classical methods. 
Of course, if the multiple factor process is stopped before the resi- 
duals are actually zero—if the matrix equation R = AA’ + U? is only 
approximately satisfied—the application of this method will give but 
approximations to the actual multiple correlations and the degree of 
the approximation will depend upon the smallness of the residuals. 

Suppose for example a factor analysis of three variables resulted 
in the same factor loadings as those of variables 1, 2, 3, in Table II, 
but with second-order residuals 712. = —.03 ; 7132 = .04 ; To32 = .02. 
Then the multiple correlation R,.; , as indicated by this method, would 
be .593. The actual value of R,..;, as determined by the correlations 
Ty» = .55, 13 = —.24, r2, = —.26 is .559 so that the approximate so- 
lution has an error whose absolute value is .034. The actual average 
of the absolute values of the residuals is .030. If the residuals had 
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been 72.2 = —.003, ris2 = .004, and 723.2. = .002, with an average ab- 
solute value of .003, the resulting error in R,..; turns out to have an 
absolute value of .0034. This subject needs further investigation, but 
it appears that there may be some definite relation between the size 
of residuals and the size of the error in the coefficient. This question 
might well be added to those recently raised by Wilson and Worces- 
ter (6, pp. 137-140) on the adequacy of reproduction of correlation 
coefficients. 

6. The case of a single factor. In the case of a single factor the 
formula (14) becomes 





(14+ 3 Bit) 
R; =: 1 ae: U;? ——__*__- —».» (16) 
\ if a >= Bi, ——T Bj; 


As an illustration, we take a one-factor problem involving eight vari- 
ables which Thurstone (5, p. 148) reproduced from Holzinger. Thur- 
stone gave the factor loadings a;,. These are given in Table III, and 
the values R; are computed from them. 


TABLE III 
Multiple Correlations in the Case of a Single Factor 


(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
Var. 





° © © 9 M,, 
j eo. a at Ae Ahan ~~ Kp 8p «= 
M, 1)/( 
1 -765 586 414 1.415 6.343 1.223 506 494 .703 
2 -739 546 454 1.203 6.555 1.184 538 462 680 
3 -716 512 488 1.049 6.709 1.156 564 436 660 
4 672 451 549 821 6.937 1.118 .614 386 621 
5 634 402 598 672 7.086 1.095 655 345 587 
6 597 357 .643 555 7.203 1.077 693 307 554 
7 095 354 .646 548 7.210 1.076 .695 305 552 
8 576 ol 669 495 7.263 1.068 -714 286 535 
7.158 = M,, 


These results are but approximations to the actual multiple correla- 
tions, since the first-order residuals (5, p. 148) are not all zero. 

The multiple correlation formula may appear in different forms 
when there is but a single factor. Using u;? = 1 — h,? and 


Qj,? as h;? 


u;" u;” : 


B;,? — 
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(16) reduces to 


1 
R;=h, |*-7 (17) 
1+ > Bi? —B;; 
é=2 


and it follows at once that R; < h; so that h; serves as an upper bound 
for R;. In general, too, as the number of variables increases without 
limit the value 1 + >B;,? — B;,? increases without limit so that R; 
converges to h; . 

This convergence of R; to h; has been surmised by Roff (7, p. 6) 
who demonstrated it for the special case in which the 7;, are all equal 
and the 7, , k ~ 1, are also all equal. His result can be obtained from 

















(17) by utilizing the fact that a2, = ds: = +++ = Qn, @j:? = Toe When 
—— v; 
7413 4,0;, = rj, so that aj, = V7, and ayn = = . We have then 
Vik 
- n—1 
=rn4/ TP @—Drp’ sie 


which is the formula given by Roff (7, p. 6) if n is replaced byn +1. 

7. The case of many factors. Approximation method. An ap- 
proximation method is desired in case there are many factors, as the 
evaluation of the 7th-order determinants demands considerable work. 
The approximation method is based on the tendency, when the num- 
ber of variables is large, for the non-diagonal terms of [J + BB’] to 
be small in comparison with the diagonal terms. This situation has 
been used previously, (8, p. 175) in determining approximate values 
of the factor loadings of. additional variables and in obtaining ap- 
proximations to the roots of the characteristic equation. In the pre- 
vious case it was assumed that Sa;,a;; = 0, while now it is assumed 
that SBixBia = Mi = @ if k a l. 

This condition is approximately satisfied by many problems even 
when the number of variables is relatively small. For example, the 
entries in column (10) of Table II are small compared with the en- 
tries in columns (9) and (11). If these terms are neglected, we get 
columns (12), (13), (14), and (15) of Table IV. Column (12) was 
obtained by multiplying the entries in columns (9) and (11) of Table 
II. Column (16), the correct correlation, is inserted to facilitate com- 
parison with column (15), the approximate correlation obtained by 
neglecting non-diagonal terms. 

The use of such approximation methods simplifies the work 
greatly if there are more than three factors. The formula becomes 


M. M. M. 
R,; ae ies 2 11 ) 22 ) M.A ( rr ) 4 
. g m (7 (w= Menic (19) 
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TABLE IV 
Approximate Solution of First Problem 
(12) (18) (14) (15) (16) 
Variable Actual 
 AU+BB’] ; 
4;[1+BB'] iS [+BB) R; R; R; 
1 57.435 1.0776 407 688 639 
2 44.157 1.0416 -750 .866 814 
3 42.157 1.4681 .560 -748 -731 
4 43.838 (1.4118) .449 .670 .668 
5 51.073 1.2118 443 .666 .666 
6 54.852 1.1283 .458 677 673 
7 36.055 1.7166 -708 841 841 
A= 61.891 


In general, as the number of variables increases, (19) gives better 

11 
M 11) j( 
is \/1 — u;? = h;. Also, if periodically there is at least one B; > e, 


Mis —>1lasn— o so that R; > h;. See also (12, Section 14). 
didi ( 


approximations. Also 





> 1 and the largest possible value of R; 


In a previous study (9, p. 10), the square root of the communal- 
ity was used as an upper bound for the multiple correlation coef- 
ficient. It is possible, with the use of (19), to obtain a good estimate 
of the correlation itself. In Table V are presented the values h; used 
previously and the new values R; obtained with (19). 


TABLE V 
Values of R; and h; for 19 Variables 


Variable 1 2 3 4 5 6 7 8 9 10 11 12 138 14 15 16 17 18 19 
h; 86 .91 .92 .94 .98 .92 .93 .95 .78 .93 .93 .90 .71 .90 .58 .89 .78 .79 .91 
R; 85 .88 .88 .91 .95 .90 .91 .93 .77 .91 .90 .89 .69 .88 .56 .87 .77 .76 .89 


8. Use of non-uniqueness methods. The development of the 
above formulas demands uniqueness factors since u; must be different 
from 0. This appears to be a serious restriction at first, and yet it is 
not so serious in practice, since observed variables usually have unique 
components (5, p. 130). However, it does not appear at first that 
methods which assume no uniqueness factor (methods which use 
AA’ = Ras the basic matrix equation) should yield results which are 
adaptable to the formulas of this paper. Of course, if the attempt is 
made to distribute the variance completely among 7 common compo- 
nents, and the attempt is successful, then the u; are zero and the for- 
mulas above are not applicable. However, this attempt is frequently 
not entirely successful, particularly if the number of observed vari- 
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ables is large, and that is the case where the method is especially val- 
uable. In an unsuccessful attempt to distribute the variance among 
r factors, the residuals of the diagonal terms are appreciably larger 
than the residuals of the non-diagonal terms. Hence, the sums of ‘the 
squares of the factor loadings are appreciably less than unity when 
the sums of the products of the factor loadings do not differ appreci- 
ably from the correlations. In such a case, even though the basic 
assumption is that u; = 0, the residuals of the diagonal terms may be 
used as the u;? and the multiple correlation coefficients may be com- 
puted by the methods of this paper. The case of perfect prediction, 
where the 2; are actually as well as theoretically zero, has been dis- 
cussed by Guttman (12, section 17). 


PART II 


PARTIAL CORRELATIONS AND REGRESSION COEFFICIENTS 
FROM THE FACTORIAL MATRIX 


1. Introduction. It has been shown in Part I how the multiple 
correlation of an observed variable with other observed variables can 
be obtained from the factorial matrix. It is the purpose of this part 
to show how partial correlation coefficients and regression coefficients 
can be obtained similarly. The general plan is similar to that of Part 
I, though the matrices involved are not so simple. We use the well- 
known determinantal formulas 





Ayo 
112.3400en — ————; 1 
2°34 VA,A_ ( ) 
A ® 
Brz-se-con — _ . (2) 
A 
Bor-saeoon a = . (3) 
It has been shown, (Part I, Equation 12) that 
Aj = UW, We ++? Wi Win ++ Ww, Ai [I+ BB] (4) 


so that A, and A. can be evaluated from the factorial matrix. The 
main task then is to evaluate 4). . 

2. The evaluation of A;.. Perhaps the simplest approach is to 
carry through an evaluation of A in which, for a number of steps, the 
value of the minor is not changed. Thus if 4,.[R] represents the de- 
terminant of the correlation matrix with the first row and second col- 
umn eliminated, then 
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AalR] = 4n| R - ae: — die He 


oe £ —A’ I I 


U: A 
ane ae [by ae the frst row 


by uw, , and the n — 2 following rows and columns by uw, --+ , Un, re- 
spectively. The result can be written 


We now divide the first column of Aa| 


0 0 B, 
Ayo = UjUgls?U,? +++ Uy? A 0 I Byix | ; 


—B’, —B' 12, I 
where B, is a single-rowed matrix consisting of = ; = P = atts - ; 
U2 U2 ~ Ue U2 


: . A Qi : 
B, is a single-column matrix composed of elements a ; Bye, is a ma- 


trix of m — 2 rows and r columns in which no one of the first sub- 
scripts is 1 or 2. B’),., is the transpose of B),.,. It follows then that 


0 0 B, 
Ayo = UjUsUs? +++ Uy? A 0 LT By 
—f’, 0 I+ B) 10, B’ 19 


0 Bz 
—_ 2...u,27A , 
U,U2U3 Un |e, I+ Byr2, Byac| 


9 BB 


= UjUsUz” +++ U,? A |, I 4 BB’ (5) 


See also (12, section 5). | 
3. The formula for partial correlation. It follows, by substitu- 
tion of (4) and (5) in (1) that 


0 Bz 


Uy, Ug,” sits a A |B, I + BB’ 
Ax. 














Ti2. Bheeen = 


VAdA. vVluu,?---u, A,(1+ BB) ][u,%,*--- th? A, + BB)] 
Pe oe 
—B, 1+ BB’ 


— e 6 
VA,[ + BB) A,[! + BB) ” 








The determinant in the numerator has r + 1 rows and columns, while 
the determinants in the denominator have but r rows and columns. 
The matrix [J -++ BB’] consists of r rows and r columns having ele- 
ments 
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The determinants of the denominator are similar to the determinant 
of [I + BB’], but the variable x, is omitted in the first case and the 
variable x, in the second. The calculation is relatively simple and fol- 
lows quickly from results used in getting the multiple correlation. 

It is readily established that, more generally 


i+j+1 0 B; 
i al _By 1488 | 
VAi[I + BB] 4;{I + BB’) © eS 


The reader should compare (12, section 7). 

4, Illustration. As an illustration we take the two-factor prob- 
lem used as an example in the case of multiple correlation. We take 
out the first three variables, for ease in checking by classical meth- 
ods, and compute 7,2... The factor loadings are augmented to form the 
B’s and the values M;,; determined. We find, from Table I of Part I, 
B,, = .889712, Bi. = .163082, B., = 1.694908, B.. = .542176 , 


M,, = 5.640949, M,,. = M2, = —.087051, M.. = 2.677230, 
with 








Tijer20n = 


A, {I+ BB’] = 12.799994, A.[1-+ BB’] = 5.585452. 
| .000000 1.694908 .542618, 


4 |B, I mn * — —.889712 5.640949 —.087051/ = 4.602031 
, —.163082 —.087051 2.677230 
so that 





= V (12.799994) (5.585452) 
The value 7,..,; as determined by the classical method is 


Ain _ S016 _ 544071 . 
V4.4, .9216 





12-3 = 


Of course, as in the case of multiple correlation, it is better to use 
classical methods if as few as three variables are involved. The ad- 
vantage of the factorial method increases with the number of vari- 
ables involved and particularly so if many partial correlations are 
desired. The flexibility, too, should be noted since it is possible to 
obtain the partial correlation between two variables against any sub- 
set of the original variables. 
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5. The case of a single factor. In the case of a single factor, 
Atl+ BB] =1+D>B,, B =B,, —Bi=—Bi,, 
Aifi+ BB] =1+)> B,— Bi, 
A,{I+ BB’] =1+ 5 B*.. — B’;, ; so that 
Rie B?,, B?); 
\ OFS SSS)’ ® 
From the one-factor problem used in Table III we note that B*;, = 


1.415, B?;, = 1.203 and, if we consider the first three variables only, 
1 + B?,, =f. B?,, = 3.252 F i ot B*,, a B?,, == 3.464, and 


_ | (1.415) (1.203) 
"we * | (S252) (8.464) 


This is in error because the residuals were not exactly 0. The actual 
value of 7,..;, as found by conventional methods, is .428. The abso- 
lute value of the difference is .039 and is about the same size as the 
average absolute value of the size of the three residuals (5, p. 148), 
which is .035. 

It can be shown in general from (8) that, in the case of one fac- 
tor, the partial correlation approaches 0 as the number of variables 
increases. The old factor loadings remain the same with each new 
addition and if periodically there is added at least one factor loading 
as large as e, the denominators increase while the numerators re- 
main the same. 

The ease with which partial correlations of high order are com- 
puted, in the case of one factor, is shown in Table I where the values 





Tijerzen = (—1) t+j+1 











OF M13, Maes » Maser » T12+345 » T12+3456 » 112634567 » Tie-za56e78 » AYE exhibited. 
TABLE I 
Computation of the partial correlations 7,,. 
Variable a; ,? u;? B,,? B,,?B,,? Vio. Mi yy 1¢ M, 1) 2 M1) 15." Ty. 
nett | 
.086 .414 1.415 1.702 
.546 .454 1.208 1.702 2.203 2.415 5.820 .3820 .566 
.512 .488 1.049 1.702 .3 8.252 3.464 11.265 .151 .389 
.451 .549 .821 1.702 .34 4.073 4.285 17.458 .0975 .312 


402 .598 .672 1.702 .345 4.745 4.957 28.521 .0724 .269 
357.643 .555 «1.702 .3456 5.300 5.512 29.214 .0583 .241 
854 .646 .548 1.702 .34567 5.848 6.060 35.439 .0480 .219 
381 .669 .495 1.702 .845678 6.343 6.555 41.578 .0409 .202 


Cn moor NW 


The 1764 different partial correlation coefficients of this prob- 
lem can be computed in this way. 
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6. Alternate form of the general formula. The numerator of 
(7) can be expressed as the sum of two determinants of the rth order 
rather than one of the (7 + 1)st order. Since 


I B, a , ye. B, 
4|_ py, I+ BB | its selistnttiatiniin |B I+ BB | 


we have 


, B, am, ; a I B, 
4|_», 14 pp |=—4U+8B1+4|_» 14 Be | 


= — A[I+ BB] +4 g 14.BB + BB.) | 
= — A[I+ BB’])+ AftI+ BB+ B.B’,]. (9) 
For example, in the illustration of section 4, 
A{I + BB'] = 15.094540, while 


[7.148929 .395329 ] __ 
aU + 83+ 81 = | 189358 2.765149 Shreeeneteth 


and 
A a = 4.602030 , as indicated earlier 
—B', I+ BB : ‘ ‘ 

7. Approximation Formulas. The computation can be consid- 
erably simplified if we make an approximation by placing the non- 
diagonal elements of [J +- BB’] equal to zero. In this case the matrix 
[J + BB’] is a diagonal matrix with elements M,,. Then 


0 B- = Mii 
A | _, rape =2 Bu Bu a 9 and 


_- By. Bu Mii / Mux (16) 
Vil Miya I Miya 
As an illustration, we compute an approximation to 112.3456, for 
the two-factor problem of section 4. In this case, if we take the fac- 
tor loadings to three places we get B,,B., = 1.506, BicBor = .0885, 
M,, = 11.753 ’ M.. = 5.266 ’ Marya = 10.963 9 My) = 8.883 , Mes) 





T1203400n 





= 5.239, Meso, = 4.971, with 
_ (1.506) (5.266) + (.0885) (11.753) __ 7 


112634567 ee si Oo eran cea 


V (10.963) (8.883) (5.239) (4.971) 


A further approximation can be made, if n is large enough so 
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that M;; , Mj;),, and M;i)., are approximately equal. In this case 
cme s Bux Bu 


T12+34-0en = 
k=1 Mix 





(11) 


It can be shown that fio.34...n > 0 as nm > o if periodically some 
ai; > ¢, since By, and By, remain the same, but Mj, increases. See also 
(12, section 15). 

8. Regression coefficients. From (2) and (5) it is apparent 
that 











a UqUoUg? +++ My? A |e, 14 Be | 
Biz-s4e00n = A, “tk Un?Us? +++ Un? A, {I + BB’] 
0 B, 
“4 |e, ee BB | 
“ A, + BBY} 


In general, similarly 


0 B; 
nibs |e I i) 
Pitan = Kak) 1 oe ey _ 


Approximations can be made as in the last section. See also (12, sec- 
tion 7). 





PART III 


THE DIFFERENT MULTIPLE CORRELATIONS INVOLVING THE 
OBSERVED VARIABLES AND THE COMMON VARIABLES 

1. Introduction. In this part a study is made of the multiple 
correlation coefficient (1) when an observed variable is estimated 
from (a) common variables, (b) common variables and observed vari- 
ables, (c) observed variables; (2) when a common variable is esti- 
mated from (a) common variables, (b) common variables and ob- 
served variables, and (c) observed variables. 

In this part, we make use of the fact that the correlation between 
an observed variable x; and a common variable y; is equal to a;; . This 
well-known fact follows at once from [Part I, Equation (1)] by form- 
ing 


= ej Yi Dy yi? _ 





ne ae =e. (1) 











P. S. DWYER 227 


The matrix of the correlations of the n observed variables and the 
r common variables can be written in the form oa r lit the observed 


variables precede the common variables horizontally to the right and 
vertically downward. 

2. The multiple correlations of an observed variable. The mul- 
tiple correlations of the observed variables have been worked out pre- 
viously. Results and references are here stated. 

(a) Correlation with the common variables. 

The multiple correlation of the observed variable x; with the r 
common variables has been shown by various authors (7, p. 2) (10, 
p. 164) to be equal to h; . 

(b) Correlations with the common variables and the n — 1 ob- 

served variables. 

It is likewise well-known (7, p. 4) (10, p. 165) that the multiple 
correlation of an observed variable with the 7 common variables and 
the n — 1 observed variables is equal to h;. 

c) Correlation with the n—1 observed variables. 

It is shown in Part I that 


[yup SUE BET 
V0 AyLT-FBBT * 


It is not necessary to use all the n — 1 other variables in apply- 
ing this formula. 

3. The multiple correlations of the common variable. 

(a) With the r— 1 other common variables. 

This trivial case results in R; = 0, since A = A; = 1 in the for- 
mula (4, p. 301). 





(2) 


| A 
\ 1— - (3) 
(b) With the r—1 other common variables and the n observed 
variables. 
The value of A in the determinantal formula (3) is now 


4| a: {| =4wa 


a= 


as shown in Part I. The value of A; is not quite so easily obtained 
since it calls for the deletion of row n + i and column n +- i in the 
matrix having n + r rows and n + r columns. Now the order in 
which the y’s occur is not material so, for convenience, we assume 
that y; is the last variable y,. We then write the factorial matrix in 
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two parts. The first part, which we designate A),,, indicates the ma- 
trix composed of all the columns of the factorial matrix with the ex- 
ception of column 7. The second part of the matrix, which we desig- 
nate by A,, consists of the rth column only. In this notation the ma- 
trix of the intercorrelations of the n + 7 variables becomes 


R A) A, 
MB yeq TF 0 
A’, 0 I 


while the determinant A, can be written 


“oe a eer U2 Ay, A, 
A A’ I “ T=AlLA's I 0 |;=4 0 TI 0 
re 0 0 I a 0 I 


wes 1+ B,B', By B, 
-atwaal a I 0 |= atv af rs 
I 











—B’, 0 0 0 I 
= A[U*] (1 + B,,+ B*., +- ++ By) 
= ALO) (14 et ot +) (4) 
It follows at once, since i may be substituted for 7 , that 
EO ie ily tie’ oc. = > Bi: 
R= 1-5 Ji-rse res 


This formula is also given by Guttman (12, 9a). It is clear that 
R; < 1 and that R; > lasn — «o if periodically some B?;; is at least 
of size e. 

The formula (5) is adapted to easy computation. It is to be noted 
that the quantities B;; are quantities which are used in finding the 
multiple correlation of an observed variable with the other observed 
variables. 

(c) With n observed variables. 

We assume again, without loss of generality, that the common 
variable is the last common variable, 7. We then have 


bey Pb. eR. ae 
&, = J1—s- with 4 = 4 sii 7 [and 





A, = A[R] = A[U*]A[I + BB’ as shown in Part I. Now 


wer “—— e. U2 Ay, A, 
A A’ 1 j — A 0 I 0 = A —A’),; I 0 
Tay See ee oe ce 
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I By B, 
= A[U?] 4B I 0 |= 4[U*?]4+ B,,, Bd. 
0 0 I 


It follows that, in general, 





R= hh — AUF By Byid 
a A(T -+ BB} 
The reader may refer to Guttman (12, Section 10). 
Here A[I + B);, B’)i,] differs from A[J + BB’] only in having 
row 7 and column 7 deleted. As an approximation formula we may set 


all non-diagonal elements in row i and column i of J + BB’ = 0 and 
we get 





(6) 





¥e 1 
B= JI-TEsR, ~ 


which agrees with the result (5) of the previous case. 

4, Illustration. As an illustration we take a problem used by 
Thurstone (5, p. 128) and compute the values of the B’s to three sig- 
nificant places. These are presented in Table I. 


TABLE I 
The Computation of the B’s. 
(1) (2) (3) (4) (5) (6) (7) (8) 
Var. 
j a; Qi h;? u;? u; Bj, Bio 
1 .660 121 .450 .550 -742 .889 168 
2 .830 .266 -760 .240 490 1.694 .543 
8 —.541 .638 -700 .300 .548 —.987 1.164 
4 —.126 771 .610 390 625 —.202 1.234 
5 487 591 .540 .460 .678 .645 .872 
6 688 —.337 .520 .480 .698 .921 —.486 
7 904 -109 .830 .170 ANZ B.104 .265 


It follows at once that 
7 
1 a > B?;; = 11.753 ’ 
j=1 
7 
= Bj, Bj. = 368 , 


1+ 3 BY, = 5.266 , 


so that, if y, is correlated with y, and the seven observed variables 
we get 





ee 
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_ [10.753 _ 
Ry = 4! 7753 ~ 96 - 


In a similar fashion 








4.266 
R,, = \ £566 =,900. 
If y, is correlated with the seven observed variables we have 
‘eae =e et 
ase = .956 
t= | eee ae 
\ .363 5.266 | 
and, similarly 
_— |: aie 
a (11.753 3.63) nisin 
\ 363 5.266 


The correlation of each observed variable with the remaining 
observed variable has been exhibited in Part I. A general summary 
of results is given in Table II. 

TABLE II 
Different Multiple Correlations. 
Predicting Variables (Excluding predicted Variable) 


Predicted Common Common Variables Observed 
Variable Variables and Observed Vars. Variables 
Z, .671 .671 .639 
Xe .872 .872 814 
L. 837 837 781 
e. .781 -781 .668 
£., .735 °735 .666 
. 721 Toi .673 
2, 911 911 .841 
Ys, 0 .956 .956 
Yo 0 .900 .900 


6. Additional Cases. Cases may arise in which one does not de- 
sire to use all the observed variables nor all the common variables. It 
is possible, in general, to work these cases by means of matrix devel- 
opments similar to those of this paper. The desired matrix is set up 
and then built up so that complete correlation and factorial matrices 
result. The determinant of this matrix is then reduced to the determi- 
nant of a matrix having but 7 rows and 7 columns. 

Frisch (11, p. 15) has used “hollow” determinants in which the 
entries in the diagonals are replaced by 0. These also can be evaluated 





a tt Ge ae ie 6h—|hCO* a 


6 #2 8 tt G&G & Mw ek he lm [CO ee 
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from the factorial matrix, but the loadings are now “communality 
augmented” rather than “uniqueness augmented” loadings. For 


ee 0 AD LI oe 
aemnaalPy! Haale! 4-48 f] 


0 I 
(8) 
.. S ‘ 
= (—1)" AL] A lc : = (—1)" AH") AUI—CC’}). 
where A[H?] = h?,h2, --- h?, and the elements of C are 
_ Ui 
Ca= i,’ (9) 


7. Conclusion. It has been shown in this paper how the results 
of a factorial matrix can be used in obtaining multiple and partial 
correlations. The presentation is illustrative rather than exhaustive, 
general rather than rigorous, but it does emphasize the actual com- 
putation of multiple and partial correlations in problems involving 
many variables. Special attention is given to formulas which, with a 
relatively small amount of work, will yield good approximations to 
the desired quantities. 

From the general theoretical point of view, it seems, there is one 
point which should be emphasized and that is that the correlation ma- 
trix results from an incomplete preliminary reduction of the data 
since the correlation matrix does not permit the isolation of a set of 
numbers which accompany one variable only. To be sure the correla- 
tion matrix does provide a set of numbers for each variable, but these 
numbers are associated with other variables; in fact, there has to be 
one of these numbers for each of the other variables. However, if the 
correlation matrix is reduced to a factorial matrix, we have each 
variable indicated by a set of numbers, numbers which are, in form 
at least, independent of the other variables, and from these numbers 
we can now construct not only the correlation matrix, but also the 
classical partial and multiple correlation coefficients as well. Since 
the numbers associated with the observed variable are not associated 
with the other observed variables, the method makes feasible the im- 
mediate study of any subset of the observed variables. This flexibility 
is most important as is the fact that the representation of a variable 
by such a set of numbers gives promise of the development of very 
objective criteria for the inclusion or exclusion of a variable in pre- 
diction (10,11, 12). It appears then that the computation of a fac- 
toria] matrix is indicated as a proper step in preliminary reduction of 
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the data by the student of multiple correlation as well as by the stu- 
dent of multiple factor analysis. 
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